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In peripheral relativistic heavy ion collisions at the Large Hadron Collider, 
jet+jet and photon+jet final states can be produced when a photon from the 
virtual photon field surrounding the nucleus interacts with a parton in the 
opposite nucleus (direct production). The virtual photon may also fluctuate 
into states with multiple gluons and qq pairs (resolved production), opening 
more channels for jet photoproduction. We compare the rates for direct and 
resolved jet+jet and photon+jet production to explore the sensitivity to the 
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Q^ ■ nuclear and photon parton distribution functions. We calculate the transverse 
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momentum distributions of both partonic jets and leading hadrons produced 



O by jet fragmentation. 
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• I. INTRODUCTION 

a: 

Ph. We discuss measurements of the nuclear parton distributions via dijet (jet+jet) and 

7+jet photoproduction in ultraperipheral heavy ion collisions at the Large Hadron Collider 
(LHC). In these collisions, the accelerated ion is surrounded by a cloud of almost real 
photons of virtuality \q 2 \ < {Tic/Ra) 2 where Ra is the nuclear radius. The virtuality, less 
than (60MeV) 2 for nuclei with A > 16, can be neglected. Since photon interactions are 
long range, they can interact with partons in the opposite nucleus even when the nuclei 
themselves do not interpenetrate. Because the photon energies are less than those of the 
nucleons, these photonuclear interactions have a smaller average center of mass energy than 
hadronic parton-parton collisions. However, even though the energy is smaller, the coherent 
photon beams have a flux proportional to the square of the nuclear charge, Z, enhancing 
the rates relative to those of photoproduction in pp collisions at the same energy. Thus 
photoproduction rates in heavy ion collisions can be high. 

Photoproduction occurs two ways in heavy ion collisions: "direct" and "resolved" pro- 
duction. "Direct" photoproduction occurs when a photon emitted from one nucleus interacts 
with a parton from the other nucleus. The photon can also fluctuate into states with mul- 
tiple qq pairs and gluons, i.e. \n(qq)m(g)) . One of these photon components can interact 
with a quark or gluon from the target nucleus ("resolved" production) [1]. The photon 
components are described by parton densities similar to those used for protons except that 
no useful momentum sum rule applies to the photon [2]. The quark and gluon constituents 
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of the photon open up more channels for jet photoproduction and could, in principle, lead 
to larger rates for resolved production in certain regions of phase space. 

In central collisions at RHIC, leading particles in jets are easier to detect above the high 
charged particle multiplicity background than the jets themselves since these high transverse 
momentum, pr, particles can be tracked through the detector [3,4]. In peripheral collisions, 
especially at LHC energies, jets should be easier to isolate and may be observed directly 
using standard high energy jet detection techniques. We thus discuss the pr distributions of 
both jets and leading particles. We work at leading order, LO, to avoid any ambiguities such 
as jet reconstruction and cone size. The ratio of the next-to-leading order, NLO, to LO jet 
cross sections appears to be relatively constant at high p T [5] . We discuss the fragmentation 
of jets and present the transverse momentum distributions of charged pions, charged kaons 
and protons/antiprotons. 



II. DUET PRODUCTION 

A. Direct dijet production 

The first hadronic reaction we study is / ~f(k) + N(P 2 ) — * jet(pi) + jet(p 2 ) + X where k and 
P 2 are the photon and nucleon four-momenta. The diagrams for direct and resolved dijet 
hadroproduction are given in Fig. 1. The two parton- level contributions to the jet yield in 
direct photoproduction are j(k) +g(x 2 P 2 ) — > q(pi) + q{P2) and r y(k) + q(x 2 P2) — > g(Pi) + q(P2) 
where also q—*q. The produced partons are massless, requiring a minimum p T to keep the 
cross section finite. At LO, the jet yield is equivalent to the high p T parton yield. The jet 
Pt distribution is modified for photoproduction from e.g. Refs. [6-8], 
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where x 2 is the fraction of the initial hadron momentum carried by the interacting parton 
and n is the momentum scale of the interaction. Recall that the extra factor of two on 
the right-hand side of Eq. (1) arises because both nuclei can serve as photon sources in AA 
collisions. The partonic cross sections are 
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The first is the photon-gluon fusion cross section, the only contribution to massive QQ 
photoproduction [9], while the second is the QCD Compton process. At LO, the partonic 
cross section is proportional to aa s (fi 2 )eQ, where a s ([i 2 ) is the strong coupling constant to 
one loop, a = e 2 /hc is the electromagnetic coupling constant, and cq is the light quark 
charge, e u = e c = 2/3 and = e s = — 1/3. The partonic invariants, s, t, and u, are defined 
as s = (k + x 2 P 2 ) 2 , t = (k- pi) 2 = (x 2 P 2 - p 2 ) 2 , and u = (x 2 P 2 - px) 2 — (k — p 2 ) 2 . In this 
case, s = 4k r y L x 2 m p where 7^ is the Lorentz boost of a single beam and m p is the proton 
mass. Since k can be a continuum of energies up to -Ebeam = lL m P , we define x 1 = k/P±, 
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analogous to the parton momentum fraction in the nucleon where Pi is the nucleon four 
momentum. For a detected parton in a nucleon-nucleon collision, the hadronic invariants 
are then S NN = (P 1 + P 2 ) 2 , T — (P 2 — Pl )\ and U = (P X - Pl ) 2 . 

The produced parton rapidities are y\ and y 2 . The jet parton rapidity is related to 
the invariant T by T = -^/~S^p T e~ yi . At LO, x Y = (pt/VSnn)^ 1 + e y ' 2 ) and x 2 = 
{PtI V 'Snn){.z~ Vi + e~ m ). We calculate x\ and x 2 as in an NN collision and then determine 
the flux in the lab frame for k = Xi'-fiirip, equivalent to the center of mass frame in a 
collider. The photon flux is exponentially suppressed for k > ^lTic/Ra, corresponding to 
a momentum fraction x\ > hc/m p RA- The maximum <yN center of mass energy, ^JS^n, is 
much lower than the hadronic ^Snn- The equivalent hadronic invariants can be defined 
for photon four momentum k as S 1 n — (k + P 2 ) 2 , T^n = (P2 — Pi) 2 , and U^n — (k — p\) 2 
[10]. The partonic and equivalent hadronic invariants for fixed k are related by s = x 2 S 1 n, 
t = UyN, and u = x 2 T^n. 

The photon flux is given by the Weizsacker- Williams method. The flux from a charge Z 
nucleus at a distance r is 



d 3 N y _ Z 2 aw 2 
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K 2 (w) + ^K 2 (w) 
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where w = kr/^L and Kq(w) and K\(w) are modified Bessel functions. The photon flux 
decreases exponentially above a cutoff energy determined by the nuclear size. In the lab 
frame, the cutoff is k max ~ ^ific/RA- In the rest frame of the target nucleus, the cutoff 
is boosted to E max = (2^ 2 — l)hc/ Ra- Table I shows, for AA collisions at the LHC, the 
nucleon-nucleon center of mass energies, y/S^N, the beam energies, -Ebeam, Lorentz factors, 
7l, maximum photon energies in the lab and rest frames, /c max and E max respectively, as 



well as the corresponding maximum jp —> qq center of mass energy, ^ S 1 n = y2E max m p , 
for single photon interactions with protons. We have also included the AA luminosities in 
Table I to aid the calculation of rates. 

The total photon flux striking the target nucleus is the integral of Eq. (4) over the 
transverse area of the target for all impact parameters subject to the constraint that the 
two nuclei do not interact hadronically [11]. This must be calculated numerically. However, 
an analytic approximation for AA collisions is given by the photon flux integrated over 
distances r > 2Ra, 
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where w R = 2kRA/lL- The difference between the numerical calculation and the analytic 
results is typically less than 15% except for photon energies near the cutoff. We use the 
more accurate numerical calculations. 

The nuclear parton densities F^(x, fi 2 , b, z) can be factorized into x and fi 2 independent 
nuclear density distributions, position and nuclear-number independent nucleon parton den- 
sities, and a shadowing function S A (x, fi 2 , b, z) that describes the modification of the nuclear 
parton distributions in position and momentum space. Then [6,12-15] 

F^x, /i 2 , b, z) = p A 0, z)S A (x, /i 2 , b, z)f l N (x, /i 2 ) (6) 
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where fi*{x, p 2 ) is the parton density in the nucleon. We use the MRST LO parton dis- 
tributions [16] at p 2 = p\. In the absence of nuclear modifications, S l = 1. The nuclear 
density distribution, p A (b,z), is a Woods-Saxon with parameters determined from electron 
scattering data [17]. Although most models of shadowing predict a dependence on the 
parton position in the nucleus, we neglect any impact parameter dependence here so that 
S' A (x,p 2 ,b,z) -> S A (x,p 2 ). 

The nuclear modification, collectively referred to as shadowing here, encompasses three x 
regions: low x shadowing where S l A < 1, x < 0.01; antishadowing where S l A > 1, 0.01 < x < 
0.1 and the EMC region where again S l A < 1 but x > 0.2. We use two parameterizations 
of nuclear shadowing which cover extremes of low x gluon shadowing. The Eskola et al. 
parameterization, EKS98, is based on the GRV LO [18] parton densities. At the minimum 
scale, yUo, valence quark shadowing is identical for u and d quarks. Likewise, the shadowing 
of u and d quarks are identical at pq. Although the light quark shadowing ratios are not 
constrained to be equal at higher scales, the differences between them are small. Shadowing 
of the heavier flavor sea, s and higher, is calculated separately at p . The shadowing ratios 
for each parton type are evolved to LO for 1.5 < p < 100 GeV and are valid for x > 10~ 6 
[19,20]. Interpolation in nuclear mass number allows results to be obtained for any input 
A. The parameterization by Frankfurt, Guzey and Strikman, FGS, combines Gribov theory 
with hard diffraction [21]. It is based on the CTEQ5M [22] parton densities and evolves 
each parton species separately to NLO for 2 < /i < 100 GeV. Although the given x range is 
10 -5 < x < 0.95, the sea quark and gluon ratios are unity for x > 0.2 and p . The EKS98 
valence quark shadowing ratios are used as input at po since Gribov theory does not predict 
valence shadowing. The parameterization is available for four different values of A: 16, 40, 
110 and 206. 

Figure 2 compares the two parameterizations for A ps 200 and p = 10, 100 and 400 
GeV. Some care must be taken when applying these parameterizations to high px since the 
upper limit of their fit range is 100 GeV. The valence ratios, identical for EKS98 and FGS 
at po, are somewhat different at high scales because the evolution is different, LO for EKS98 
and NLO for FGS. The sea quark ratios are rather dissimilar. The EKS98 ratios have no 
antishadowing but, instead, show a peak at S A 1 when x ~ 1, decreasing smoothly at 
larger x, the EMC region. The FGS parameterization has an antishadowing peak at x ps 0.07 
after which S A ~ 1 for x > 0.2, as at p . Thus quark shadowing in the EMC x region will 
be stronger for EKS98. At the lowest x and p shown, the FGS gluon ratio shows stronger 
shadowing and larger antishadowing than EKS98. By p > 100 GeV, the results are rather 
similar. It is notable that at p the FGS gluon ratio is pa 1 for x > 0.2 and decreases 
below 1 at higher scales, unlike the sea quark ratios, exhibiting an EMC-like effect. At 
high pt, therefore, the EKS98 and FGS gluon ratios should not be significantly different. 
As we have already pointed out, although the upper limit on the range of both shadowing 
parameterizations is 100 GeV, the parameterizations do not strongly deviate from the trends 
at lower scales. Thus the results can be taken as indicative of the expected behavior. 

There are a few photon parton distributions available [23-29]. The data [30,31] cannot 
definitively rule out any of these parton densities. As expected, F^(x, p 2 ) = F?(x, p 2 ) flavor 
by flavor because there are no "valence" quarks in the photon. The gluon distribution in 
the photon is less well known. We use the GRV-G LO set [23,24]. At p T > 10 GeV, the 
difference in results due to the photon parton densities is minimal. 
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The direct jet photoproduction p? distributions are given in Fig. 3 for AA interactions 
at the LHC For Pb+Pb collisions at y/S^N — 5.5 TeV, we show the px distributions of 
the produced quarks, antiquarks and gluons separately, as well as their sum. For Ar+Ar 
collisions at y/S NN = 6.3 TeV and O+O collisions at y/S NN = 7 TeV, we show only the 
total px distributions. All the results are shown in the rapidity interval |yi| < 1. 

Quarks and antiquarks are produced in greatest abundance, with only a small difference 
at high pt- Photon-gluon fusion alone produces equal numbers of quarks and antiquarks. 
The quark excess arises from the QCD Compton diagram which also produces the small 
final-state gluon contribution. The j(q + q) contribution grows with p T since the valence 
quark distributions eventually dominate production, as shown in Fig. 3(b) where the yg 
contribution is compared to the total. At low px, the contribution is ~ 90% of the 
total, dropping to 10 — 30% at pr ~ 400 GeV. At the large values of x needed for high 
Pt jet production, v > f®. Thus the QCD Compton process eventually dominates dijet 
production, albeit in a region of very low statistics. The •yg contribution is larger for the 
lighter nuclei since the higher energies reduce the average x values. 

The direct dijet photoproduction cross section is significantly lower than the hadropro- 
duction cross section. Some of this reduction is due to the different couplings. The pho- 
toproduction rate is reduced by a factor of ae^/a s ~ 100. There are also fewer diagrams 
available for jet photoproduction relative to all 2 — > 2 scatterings in hadroproduction. In 
addition, the gg — > gg hadroproduction process, with its high parton luminosity, has no 
direct photoproduction equivalent. 

Since the typical scales of jet production are large, the effects of shadowing, reflected in 
R(pt) = (d(r[S' l A ] I dpx) I \da\S\ = l]/dpr), are rather small, see Fig. 3(c) and (d), because 
the average x is high. The differences between the two shadowing parameterizations are 
on the few percent level. At low p T , the produced quarks and antiquarks are mainly from 
gluons. The produced gluons only come from quarks. The peak for the produced quarks and 
antiquarks in R(pr) between 50 < < 100 GeV is due to gluon antishadowing. The total 
R{Pt) for all produced partons in Pb+Pb collisions is dominated by the 75 contribution. 
The maximum value of S\ in the antishadowing region is ~ 1.07 for EKS98 and ~ 1.1 for 
FGS, reflecting the high fx 2 behavior of the shadowing parameterizations. 

The EKS98 ratios for the produced quarks and antiquarks in Fig. 3(c) follow R(j>t) for 
the total rather closely over all pt- The quark and antiquark ratios are slightly above the 
total at low pt due to the small 75 contribution. They continue to follow the total at high 
Pt since all the EKS98 ratios exhibit similar behavior at large x. The produced gluon ratio 
follows the quark ratios for p T > 200 GeV. The large p T contribution arises from the valence 
quarks. Some antishadowing is observed at low p T due to the valence quark contribution. 
The total ratios for the lighter ions are closer to unity for all pt due to their smaller A. 

The results are similar for FGS, shown in Fig. 3(d), but there are some subtle differences. 
The ratio R{pt) for produced gluons, arising from the yq contribution, exhibits a larger 
antishadowing effect on R{j>t) because S\ is higher for this parameterization, see Fig. 2. 
The FGS antiquark shadowing ratio goes to unity for x > 0.2, flattening R(j>t) for antiquarks 
(dotted curve) due to the contribution from 75 — > gq. The FGS valence quark ratio is taken 
from EKS98, resulting in the similarity of R(pt) i n Figs. 3(c) and (d) at high p?. 

Recall that some care must be taken when applying these parameterizations to high pt 
since the upper limit of their fit range is 100 GeV. While no extraordinary effects are seen 
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in their behavior beyond this scale, the results should be taken as indicative only. Finally, 
we remark that we have only considered the range \y±\ < 1. Including contributions from 
all rapidities would increase the effect of shadowing since smaller x values can be reached 
when | yi | is large. 

In Fig. 4, we present the rapidity distributions with two different values of the minimum 
Pt, Pt > 10 GeV on the left-hand side and 100 GeV on the right-hand side. The results, 
given by the solid curves, are shown without nuclear shadowing effects. Note that, in this 
case, the photon comes from the left. There is a symmetric case where the photon comes 
from the right, the factor of two on the p T distribution in Eq. (I). In this case, the yi 
distribution is reflected around y x = 0. With the 10 GeV cut, the distributions are rather 
broad at negative y± where the photon has small momentum and, hence, large flux. At 
large y\ > 0, corresponding to small x for the nucleon momentum fractions and high photon 
momentum, the distributions fall rapidly since at high photon momenta, the photon flux is 
cut off as k — > k mSLX . The distributions with the 100 GeV cutoff are narrower because the 
edge of phase space is reached at lower values of y 1 . The rapidity distributions are broader 
in general for the lighter systems due to the higher y/S^N. 

Figure 5 gives the ratio R(yi) = {do-[S)^ / dyi) / \da\S\ = \}/dyi) for the two px cuts. The 
ratios reflect the direction of the photon, showing an antishadowing peak at y\ ~ —3, an 
EMC region at y\ < —4 and a shadowing region for yi > —0.5 for p T > 10 GeV, the left-hand 
side of Fig. 5. The shadowing effect is not large, (20 — 25)% at y\ ~ 4 for Pb+Pb collisions 
and decreasing with A. The antishadowing peak is higher for FGS while its shadowing effect 
is larger at positive y lt as also noted in the py-dependent ratios in Fig. 3. A comparison of 
the average effect around \yi\ < 1 with the pt ratios shown in Fig. 3, are in good agreement. 
Even though x% is smaller for the lighter systems, the shadowing effect is also reduced. Since 
shadowing also decreases with /z 2 , the effect is even smaller for p T > 100 GeV, only ~ 5% 
at y± > 0, as shown on the right-hand side of Fig. 5. Here the rise at y\ < —3.5 is the Fermi 
motion as x 2 — > 1. At yi > —1.2, the antishadowing region is reached. 

The total dijet photoproduction cross sections without shadowing and the corresponding 
rates in a one month, 10 6 s, LHC run are shown in Table II. The rates are based on the 
luminosities of Table I. All results are given in the rapidity interval \y±\ < 1. Extended 
rapidity coverage, corresponding to e.g. \yi\ < 2.4 for the CMS barrel and endcap systems, 
could increase the rates by a factor of ~ 2. (The increase in rate with rapidity acceptance 
is not linear in because the rapidity distributions are asymmetric around y± = and 
increasing pt narrows the rapidity distribution. The effect of changing the y\ cut is closer 
to linear at low pt and larger at high pt because the peak is at y± < —1 for large pr, as seen 
in the y\ distributions on the right-hand side of Fig. 4.) Note that by pt ~ 100 GeV, only 
a few events are expected per month. However, while high pt jets are rare, they should be 
easily observable in a clean environment like photoproduction. 

There is a difference of ~ 500 in the Pb+Pb and O+O cross sections at pt ~ 11 GeV, 
decreasing to less than a factor of four at m 400 GeV. The difference decreases with p T due 
to the larger phase space available at high p T for the higher \/Snn systems. Note that the 
rates are nearly the same for all systems because even though A is larger for Pb+Pb, higher 
luminosities and higher v^Saw compensates for the lower A in lighter systems. 

We point out that the jet hadroproduction rates are much higher because they are 
multiplied by ~ A 2 for hard processes which increase with the number of binary collisions 
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in A A collisions. (The relation is not exact due to shadowing.) The integration over all 
impact parameters leads to ~ A 2 scaling while there is only a factor of A in the dijet 
photoproduction rate since the photon flux is already integrated over impact parameter for 
b > 2R A . This, combined with the lower effective energy and fewer channels, considerably 
reduces the photoproduction rates relative to hadropro duct ion. 

B. Direct leading hadron production 

We now turn to a description of the final-state particle production in the hadronization 
of jets. The particle with the highest pt is called the "leading" particle. The corresponding 
leading particle pt distribution is [32] 

aa ~fA-^h} 

dpT(Pb 



where the X on the left-hand side includes all final-state hadrons in addition to h but X' 
on the right-hand side denotes the unobserved final-state parton. The subprocess cross 
sections, da/dt, are related to s 2 dajdtdu in Eq. (1) through the momentum-conserving 
delta function 5(s + t + u) and division by s 2 . The integrals over rapidity have been replaced 
by an integral over center-of-mass scattering angle, 6* min < 9 cm < 6* max , corresponding to 
a given rapidity cut. Here 9 min = and # max = n covers the full rapidity range while 
#mm = tt/4 and # max = 37r/4 roughly corresponds to \yi\ < 1. The fraction of the final 
hadron momentum relative to that of the produced parton, z c , appears in the fragmentation 
function, D h /i(z c , fi 2 ), the probability to produce hadron h from parton I. The fragmentation 
functions are assumed to be universal, independent of the initial state. 

The produced partons are fragmented into charged pions, charged kaons and pro- 
tons/antiprotons using LO fragmentation functions fit to e + e~ data [33]. The final-state 
hadrons are assumed to be produced pairwise so that n = (n + + n~)/2, K = (K + + K~)/2, 
and p = (p + p)/2. The equality of p and p production obviously does not describe low 
energy hadroproduction well. As energy increases, this approximation may become more 
reasonable. The produced hadrons follow the parent parton direction. We have used the LO 
KKP fragmentation functions [33]. The KKP scale evolution is modeled using e + e~ data 
at several different energies and compared to pp, ^p and 77 data. After some slight scale 
modifications [34] all the h~ data could be fit. However, there are significant uncertainties 
in fragmentation when the leading hadron takes most of the parton momentum [35], as is 
the case here. 

We assume the same scale in the parton densities and the fragmentation functions, 
fi 2 = p\. A larger scale, p^/z 2 , has sometimes been used in the parton densities. At high 
Pt, where z c is large, the difference in the results for the two scales is small. We have not 
included any intrinsic transverse momentum broadening in our calculations [36,37]. This 
"Cronin" effect can be important when px is small but becomes negligible for transverse 
momenta larger than a few GeV. 
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The largest contribution to the total final-state charged particle production is charged 
pions, followed by kaons. The proton contribution is smallest even though at RHIC the p/ir 
ratio approaches unity in central and 0.4 in peripheral Au+Au collisions [38] for p? ~ 2 
GeV. This discrepancy is due to the poor knowledge of the fragmentation functions at large 
z c , particularly for baryons. Hopefully by the time the LHC begins operation, updated 
fragmentation functions incorporating pp data from RHIC will be available, allowing better 
estimates of leading particle production at higher z c . 

The corresponding hadron distributions from direct jet photoproduction are shown in 
Fig. 6(a) for A A collisions. Note that the leading hadron cross sections are lower than the 
partonic jet cross sections, compare Figs. 3(a) and 6(a). Several factors can account for 
this. The maximum S^n is a factor of five or more less than \/Snn for AA collisions. 
The reduced number of processes available for direct dijet photoproduction is a significant 
contribution to the decrease. Note also that the px distribution is steeper for leading hadrons 
than for the jets, as may be expected since the effective px of the hadron is higher than than 
of the produced parton. 

The average z c for direct photoproduction of high p T particles is « 0.4 for particles with 
p T 10 GeV, increasing to (z c ) > 0.45 — 0.55 for p T > 100 GeV. The lower z c values 
correspond to lighter ion collisions. In this z c region, the fragmentation functions are not 
very well known. As pointed out in Ref. [35], a small change in the fragmentation function 
fits can produce significant changes at large z c . This region is not well constrained by the 
e + e _ data used in the fits. 

The effect of fragmentation on the production channels is shown in Fig. 6(b) where 
we present the fraction of leading hadron production from the 75 channel for all charged 
hadrons. The ratios are rather similar to those of the partonic jets although the 75 fraction 
is somewhat smaller due to the larger average x of hadron production with respect to jets, 
as we discuss later. 

The shadowing ratios for charged pions produced in Pb+Pb collisions by quarks, anti- 
quarks, gluons and the total from all partons, are shown for the EKS98 and FGS parameter- 
izations in Fig. 6(c) and (d). The ratios for pion production in Ar+Ar and O+O collisions 
are also shown. The high p T FGS antiquark ratios flatten out relative to the EKS98 ratio 
because the jq channel dominates gluon production at high p T . The flattening begins at 
lower pt here because the x for hadron production is higher than that for the jets. The ratio 
of pions arising from produced gluons follows the valence ratio, as expected. The ratios 
decrease with increasing p T due to the EMC effect for x > 0.2 when p T > 100 GeV. 

We now discuss the relative values of the nucleon momentum fraction, x for parton 
and hadron production. On the left-hand side of Fig. 7, we compare the average x values 
for produced gluon jets (upper plot) and for pions produced by these gluons (lower plot). 
We have chosen to compute the results for produced gluons alone to better compare with 
resolved jet photoproduction, discussed next. Since we are interested in produced gluons, 
we only consider the QCD Compton contribution, r yq — > gq. This channel is biased toward 
larger momentum fractions than 75 —> qq since the gluon distribution is largest at small x 
while the valence quark distribution in the proton is peaked at x ~ 0.2. The average x for 
a gluon jet is ~ 0.005 - 0.008 at p T « 10 GeV, increasing to ~ 0.03 - 0.05 at 50 GeV. The 
smallest x is from the highest energy O+O collisions. The average x increases with p?, to 
~ 0.25 — 0.4 at p T ~ 400 GeV. When final state pions are considered, in the lower left-hand 
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plot, at low pt, the average x is larger than for gluon jets. At pt ~ 10 GeV, (x) 0.02 — 0.03 
while at 50 GeV, (x) « 0.09 — 0.12. At high p^, however, the average x becomes similar for 
jet and hadron production as (z c ) increases to 0.6 — 0.7 at p T ~ 400 GeV. 



C. Resolved dijet production 

We now turn to resolved production. The hadronic reaction, 7iV — > jet + jet + X, is 
unchanged, but in this case, prior to the interaction with the nucleon, the photon splits into 
a color singlet state of qq pairs and gluons. On the parton level, the resolved LO reactions are 
e.g. g(xk) + g(x2P2) — > g{Pi) + g{P2) where x is the fraction of the photon momentum carried 
by the parton. The LO processes for resolved photoproduction, illustrated in Fig. l(c)-(h), 
are the same as those for LO 2—^2 hadropro duct ion except that one parton source is a 
photon rather than a nucleon. The resolved jet photoproduction cross section for partons of 
flavor / in the subprocess ij — > kl in AB collisions is, modified from Refs. [6-8], 
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ij= 
(kl) 

1 



X 



1 + 5, 



fci 



' ^da ij ^ kl ^ r ^da ij ^ kl ^ -/ 
^ (t,M)+(5flS —j—(u,t) 
dtdu dtdu 



(8) 



where s = (x/c + X2-P2) 2 , i = — pi) 2 , and -u = (x 2 P2 — Pi) 2 - The 2 —» 2 minijet 
subprocess cross sections, da/dt, given in the review article by Owens [39], are related to 
da /dtdu through the momentum-conserving delta function 8(s + t + u). The sum over initial 
states includes all combinations of two parton species with three flavors while the final state 
includes all pairs without a mutual exchange and four flavors (including charm). The factor 
1/(1 + Ski) accounts for identical particles in the final state. 

The resolved jet results, shown in Fig. 8, are independent of the photon parton densities 
for pt > 10 GeV. Along with the total quark, antiquark and gluon cross sections in Pb+Pb 
collisions, we also show the individual partonic contributions to the jet pr distributions. 
The produced gluon contribution dominates for p T < 25 GeV but, by 50 GeV, quark and 
antiquark production becomes larger due to the increase of the qg — > qg channel relative 
to the gg — > gg channel. We also show the total px distributions for Ar+Ar and O+O 
collisions. For lighter nuclei, the crossover between gluon and quark/ antiquark dominance 
occurs at higher p T due to the higher collision energy. 

Note that the resolved dijet photoproduction contribution is larger than the direct by a 
factor of two to three for p T < 50 GeV, despite the lower effective center of mass energy of 
resolved production. Resolved production opens up many more channels through the parton 
components of the photon. Indeed, now all the 2-^2 channels for LO jet hadroproduc- 
tion are available. In addition, the quark and antiquark distributions in the photon are the 
same. These distributions are large at high momentum fractions, higher than the quark and 
antiquark distributions in the proton. Thus the quark and antiquark channels are enhanced 
relative to hadroproduction. The largest difference between the quark and antiquark pro- 
duction rates is due to the difference between the valence and sea quark distributions in the 
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nucleus. Where the valence and sea quark contributions are similar, as for \yi\ < 1, the 
difference is rather small. If all rapidities were included, the relative quark and antiquark 
rates could differ more. 

The direct and resolved rapidity distributions are compared in Fig. 4 for the two pt cuts, 
10 and 100 GeV. While the \y±\ < 1 resolved contribution is a factor of two to three larger 
than the direct at pt < 50 GeV, a comparison of the y\ distributions over all rapidities shows 
that the resolved contribution can be considerably larger, a factor of ~ 5 — 10 at y± < —3.5 
for p T > 10 GeV. At pt > 100 GeV, the resolved contribution is still equivalent to or slightly 
larger than the direct at y x < —3 but drops below at larger rapidities. Thus, going to higher 
Pt can separate direct from resolved production, especially at forward rapidities. Recall that 
the produced gluons dominate resolved production at pt < 25 GeV while they are only a 
small contribution to direct production. The largest gluon production channels are typically 
99 99 an d gq — > gq- As y\ becomes large and negative, the photon x decreases while 
x 2 of the nucleon increases, leading to the dominance of the gq channel. The photon gluon 
distribution is largest as x decreases. The valence quark distribution of the proton is also 
important at high pt, causing the resolved to direct ratio to flatten for yi > —2.5 when 
p T > 100 GeV. 

In Fig. 5, we compare the direct and resolved shadowing ratios, R(yi)- A smaller shadow- 
ing effect is observed for the resolved component due to the larger x 2 for resolved production. 
The difference in the direct and resolved shadowing ratios is reduced for larger p T - 

To measure the nuclear parton densities most directly, it is preferable for direct produc- 
tion to be dominant. However, Fig. 9 shows that a pt cut is not very effective for dijet 
production, as previously mentioned, even at forward rapidity. The resolved to direct pro- 
duction ratios are all larger than unity for p T > 10 GeV, even for large, positive y 1 . While 
the ratio is less than 1 for y x > —2.5 and p T > 100 GeV, it is only ~ 0.5 for Pb+Pb, 
increasing to 0.8 for O+O. 

Thus, although clean separation is possible at pt > 100 GeV, precision parton density 
measurements are not possible at these values of pt due to the low rate. Other means 
of separation must then be found. Resolved processes will not be as clean as direct in 
the direction of the photon due to the breakup of the partonic state of the photon. The 
multiplicity in the photon fragmentation region will be higher than in direct production 
where the nucleus should remain intact. A cut on multiplicity in the photon direction may 
help separate the two so that, although there should be a rapidity gap for both direct and 
resolved photoproduction, the gap may be less prominent for resolved production. 

The total resolved dijet photoproduction cross sections without shadowing and the cor- 
responding monthly (10 6 s) LHC rates are given in Table III. At low pt, the cross sections 
and rates are a factor of 2 — 3 higher than for direct dijet production. The largest increase is 
for the lightest nuclei since the lowest x values are probed. However, with increasing pt, the 
phase space is reduced. The average photon momentum is increased and, at large photon 
momentum, the flux drops faster. The average momentum fractions probed by the nuclear 
parton densities grows large and only valence quarks contribute. The lower effective energy 
of resolved relative to direct photoproduction reduces the high pt phase space for resolved 
production. Thus, at the highest pt, the rate is reduced relative to direct by a factor of 4 — 9 
with the smallest decrease for the lightest system due to the higher effective y/S^N- Since 
resolved production has a narrower rapidity distribution than direct production, increasing 
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the rapidity coverage would not increase the rate as much as for direct photoproduction. 

We show the individual partonic shadowing ratios for Pb+Pb collisions with the EKS98 
parameterization in Fig. 8(b). The quark and antiquark shadowing ratios are very similar 
although the quark ratio becomes larger for higher p T (higher x) values of x due to the 
valence distribution. Now the gluon ratio shows larger antishadowing since gluon production 
is dominantly through the gg and qg channels rather than / yq in direct production, compare 
Fig. 3. The FGS parameterization gives similar results, Fig. 8(c). However, since the small 
FGS gluon antishadowing is stronger, R(pr) is larger for p T < 150 GeV. 

D. Resolved leading particle production 

We now turn to leading particles from resolved jet photoproduction. The leading particle 
Pt distributions are 

dv^A^hx = 4pT f dz f e — d6 cm r°° dk d 3 N y dx dx 2 

dp T d 2 b J ie min sin 6 cm ife min k dkd 2 b Jk min /k x Vx 2min x 2 

x E ^)Ff{x 2 , n\ b, z) + Fj(x, n 2 )F t A (x 2 , /i 2 , b, z)} 

ij= 
(kl) 

1 

x l + S k i 

The subprocess cross sections, da/di, are related to s 2 da/didu in Eq. (8) through the 
momentum-conserving delta function S(s + t + u) and division by s 2 . The drop in rate 
between jets and high p? hadrons is similar to that in direct photoproduction, as can be 
seen by comparison of Figs. 8 and 10 relative to Figs. 3 and 6. Now that gluon fragmen- 
tation makes a larger contribution, the relative pion contribution is larger than in direct 
photoproduction while the relative proton contribution is significantly reduced. The smaller 
effective center of mass energy for resolved photoproduction lowers the phase space available 
for fragmentation. Baryon production is then reduced compared to light mesons. 

The reduction in phase space for leading hadrons relative to fast partons can be seen 
in the comparison of the average x values for resolved photoproduction of jets and leading 
hadrons, shown on the right-hand side of Fig. 7 for gluons and pions from gluons respectively. 
At p T ~ 10 GeV, the average x of the gluon jet is 0.03 — 0.04, increasing to 0.16 — 0.24 
at px ~ 200 GeV, higher than for direct photoproduction, as expected. The x values for 
hadron production are larger still, ~ 0.06 at low pt while x ~ 0.23 — 0.33 at pt ~ 200 GeV. 

The shadowing ratios in Fig. 10 also reflect the increasing x. Now the antishadowing 
peak is shifted to pt ~ 30 GeV since the average x values are in the EMC region, even at 
low p T - The values of R(j>t) at high p T are somewhat lower than for direct production due 
to the higher x. The average z c of the fragmentation functions is also somewhat larger for 
resolved production, 0.7 — 0.8 at pt ~ 400 GeV. 

Since the resolved jet cross section is larger than the direct at low pt, it is more difficult 
to make clean measurements of the nuclear gluon distribution unless the two contributions 
can be separated by other methods. Instead, the large valence quark contribution at high 
Pt suggests that jet photoproduction can probe the nuclear valence quark distributions at 
larger /z 2 than previously possible. At pt > 100 GeV, more than half of direct jet production 
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°>fk- 



dt 



,1 m t d(jii 
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dt 
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is through the channel. Unfortunately, the rates are low here, making high precision 
measurements unlikely. However, the events should be very clean. 



III. 7+ JET PRODUCTION 
A. Direct 7+jet production 

A clean method of determining the quark distribution in the nucleus at lower pt is the 
process where a jet is produced opposite a photon in the final state, Compton scattering in 
direct production. The cross sections are reduced relative to the jet+jet process since the 
coupling is a 2 eg in the coupling rather than aa s eg, as in dijet production. In addition, the 
quark distributions are lower than the gluon, also reducing the rate. The diagrams for 7+jet 
production are shown in Fig. 11. 

We now discuss the jet and leading particle distributions for direct and resolved 7+jet 
photoproduction. The hadronic process is 'y(k) + N(P2) — > 7(^1) + jet(p2) + X. The only 
partonic contribution to the 7+jet yield in direct photoproduction is ^{k) + g(x2-P2) — ► 
7(Pi) + q(P2) (or q — > q) where the produced quark is massless. We now have 

;2^dir 



jd 

=q,q 

where the partonic cross section for the Compton process is 



d (J^wri — >~va 2 



,2" u 7<2— H I _ ^^„,2„4 

dtdu 



s — — = --7ra e Q 



s 2 + u 21 



su 



5(s + t + u). (11) 



The extra factor of two on the right-hand side of Eq. (10) again arises because both nuclei 
can serve as photon sources in AA collisions. The kinematics are the same as in jet+jet 
photoproduction, described in the previous section. 

The direct 7+jet photoproduction results are given in Fig. 12 for AA interactions at the 
LHC. We show the transverse momentum, pt, distributions for all produced quarks and 
antiquarks in Pb+Pb, Ar+Ar and O+O collisions for \yi\ < 1. 

There is a drop of nearly three orders of magnitude between the dijet cross sections in 
Fig. 3 and the 7+jet cross sections in Fig. 12. Most of this difference comes from the relative 
couplings. The rest is due to the reduced number of channels available for direct 7+jet 
production since more than half of all directly produced are gluon-initiated for p T < 100 
GeV, see Fig. 3(b). 

We have not distinguished between the quark and antiquark initiated jets. However, 
the quark-initiated jet rate will always be somewhat higher due to the valence contribution. 
When pt < 100 GeV, the quark and antiquark jet rates are very similar since x is still 
relatively low. At higher p?, the valence contribution increases so that when p? = 400 GeV, 
the quark rate is 1.5 times the antiquark rate. Since the initial kinematics are the same for 
7+jet and jet+jet final states, the average momentum fractions for 7+jet production are 
similar to those shown for the jq — > gq channel in Fig. 7. 

The shadowing ratios shown in Fig. 12(b) are dominated by valence quarks for p? > 100 
GeV. The FGS ratio is slightly higher because the EKS98 parameterization includes sea 
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quark shadowing. The effect is similar to the produced gluon ratios, at the same values of 
x in Fig. 3(c) and (d), since the final-state gluons can only come from quark and antiquark 
induced processes. 

We next present the rapidity distributions for the same two p T cuts used for dijet pho- 
toproduction in Fig. 13. Note that the rapidity distribution for p T > 10 GeV is broader at 
negative y\ than the dijet distributions in Fig. 4 because direct dijet production is dominated 
by ig — > qq at these Pt while the valence distribution entering the 7g — > 7g does not drop 
as rapidly at large x 2 as the gluon distribution. When the turnover at large negative y\ 
occurs, it is steeper than for the dijets. However, it drops even more steeply at forward yi 
because the quark distribution is smaller than the gluon at low x 2 . When p T > 100 GeV, 
the 7+jet y\ distribution is narrower than the dijets since the quark distributions drop faster 
with increasing x 2 at high pj>. 

The shadowing ratios as a function of y\ are shown in Fig. 14. They exhibit some inter- 
esting differences from their dijet counterparts in Fig. 5 because of the different production 
channels. At pr > 10 GeV, the antishadowing peak is lower at y\ ~ —2.5 and the shadowing 
is larger at yi > 0. Although this may seem counterintuitive, comparing the valence and sea 
quark shadowing ratios in Fig. 2 can explain this effect. Valence antishadowing, the same 
for EKS98 and FGS, is not as strong as that of the gluon. The sea quarks have either no 
antishadowing (EKS98) or a smaller effect than the valence ratios (FGS). Thus antishad- 
owing is reduced for direct 7+jet production. At large y±, the x 2 values, while smaller than 
those shown in Fig. 5(a) for |yi| < 1, are still moderate. Since the evolution of the gluon 
distribution is faster with fi 2 , sea quark shadowing is actually stronger than gluon shadowing 
at pt > 10 GeV and low x 2 . When pt > 100 GeV, the Fermi momentum peak is not as 
prominent because the sharp increase in the valence and sea shadowing ratios appears at 
higher x 2 than for the gluons, muting the effect, particularly for the lighter systems. 

The cross sections and rates in a one month, 10 6 s, run, shown in Table IV, are lower 
than those for hadroproduction. Direct 7+jet photoproduction proceeds through fewer 
channels than hadroproduction where the LO channels are gq — > '-fq and qq — > g'-f, the 
same diagrams for resolved 7+jet photoproduction. This, along with the lower effective 
energy and correspondingly higher x, reduces the photoproduction cross sections relative 
to hadroproduction. The lower A scaling for photoproduction also restricts the high p T 
photoproduction rate. 



B. Direct 7+hadron production 



We now turn to a description of final-state hadron production opposite a photon. The 
leading particle pr distribution is [32] 



rirr dir 

dpT(Pb 



Ap T / dz 



d6 r 



Sllll 



d 3 N 
dk 7 



dx 2 
x 2 



X 



dkd 2 b 
da^yi D h /i(z c , jj 2 ) 



(12) 



t=q,q 



dt 



Z, 



where X on the left-hand side includes the final-state gluon. On the partonic level, both the 
initial and final state partons are identical so that parton i fragments into hadron h according 
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to the fragmentation function, Dh/i(z c , fi 2 ). The subprocess cross sections, da/dt, are related 
to s 2 da jdtdu in Eq. (10) through the momentum- conserving delta function 5(s + t + u) and 
division by s 2 . Our results, shown in Fig. 15, are presented in the interval \y\\ < 1. 

The cross section for 7+hadron production are, again, several orders of magnitude lower 
than the dijet calculations shown in Fig. 6(a). At the values of z c and x important for 
dijet production, the final state is dominated by quarks and antiquarks which fragment 
more frequently into charged hadrons than do gluons. While 75 — > qq produces quarks and 
antiquarks with identical distributions, the contribution from the 'jq — > qg channel makes 
e.g. pion production by quarks and antiquarks asymmetric. We also note that for p T < 100 
GeV, 60% of the dijet final state particles are pions, p=* 33% kaons and p=* 7% protons. As 
Pt increases, the pion and proton contributions decrease slightly while the kaon fraction 
increases. In the case of 7+hadron final states, there is no initial state gluon channel. Thus 
the valence quarks dominate hadron production and the relative fraction of produced pions 
increases to 66%. The kaon and proton fractions are subsequently decreased to Pi 28% and 
6% respectively. 

The shadowing ratios, shown in Fig. 15(b) and (c) for produced pions, kaons and protons 
separately for Pb+Pb as well as the total ratios for Ar+Ar and O+O collisions, reflect the 
quark-initiated processes. We show the results for all charged hadrons here since we do 
not differentiate between quark and antiquark production. The ratios, almost identical for 
produced pions, kaons and charged hadrons, are quite different from the ratios shown for 
pion production by quarks and antiquarks in Fig. 6(c) and (d) since these pions originate 
from initial-state gluons and thus exhibit antishadowing. The results are similar to pions 
from gluon jets in Fig. 6. However, in this case the ratios are slightly higher due to the 
relative couplings. The proton ratios are lower than those for pions and kaons due to the 
nuclear isospin. The dominance of d valence quarks in nuclei reduces the proton production 
rate since d quarks are only half as effective at producing protons as u quarks in the KKP 
fragmentation scheme [33]. This lower weight in the final state reduces the effectiveness 
of proton production by the initial state, decreasing the produced proton shadowing ratios 
relative to pions and kaons. Valence quarks dominate the observed final state shadowing at 
these larger values of x, as in Fig. 7. 



C. Resolved 7+jet production 



Now we turn to resolved production of 7+jet final states. The resolved jet photopro- 
duction cross section for partons of flavor / in the subprocess ij — > k'j in AB collisions is 
modified from Eq. (8) so that now 
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The resolved diagrams are those for hadroproduction of direct photons, qg — > q*y and qq 
g~f. The 2 — * 2 minijet subprocess cross sections are [39] 



^ 2 d (Jgg 1 2 

s — t-^- = — ira s ae n 
dtdu 3 Q 

J2 



s 2 + u 2 



~2 ^ ^99 ^ 2 

s — — = -na s ae n 
dtdu 9 Q 



i 2 + u 2 ' 



S(s + t + u) (14) 



£(s + * + u). (15) 

tu J 

Note that there is no factor 1/(1 + Ski), as in Eq. (8), since there are no identical particles 
in the final state. 

The resolved 7+jet results are shown in Fig. 16 using the GRV LO photon parton den- 
sities. Along with the total partonic rates in Pb+Pb collisions, we also show the individual 
partonic contributions to the jet pr distributions in Fig. 16(a). The total yields are slightly 
higher for the resolved than the direct contribution where only one channel is open and 
the coupling is smaller. Quark and antiquark production by the qg process is dominant for 
Pt < 40 GeV but, at higher p T , gluon production dominates from the qq channel. The large 
values of x again makes the valence quark contribution dominant at higher p?. The total 
Pt distributions for Ar+Ar and O+O collisions are also shown. 

The strong antishadowing in the produced quark and antiquark ratios in Fig. 16(b) 
and (c) comes from the qg channel. The antiquark ratio is higher because the qg parton 
luminosity peaks at higher x than the qg luminosity and at lower x the gluon antishadowing 
ratio is larger. The difference between the quark and antiquark ratios increases with p? 
since the average x and thus the valence quark contribution also grow with p T . At high p T , 
the flattening of the FGS quark and antiquark ratios is due to the flattening of the gluon 
parameterization at x > 0.2. 

The final-state gluon ratio shows little antishadowing since it arises from the qq channel. 
The antishadowing in the EKS98 ratio is due to the valence quarks while the higher ratio 
for FGS reflects the fact that the antiquark ratios also show antishadowing for x < 0.2. The 
ratio for the total is essentially the average of the three contributions at low pr, where they 
are similar, while at high p T , where the qq channel dominates, the total ratio approximates 
the produced gluon ratio in both cases. 

The resolved rapidity distributions are also shown in Fig. 13 for the two pr cuts. The 
resolved distribution is not as broad at negative y\ as that of the dijet process in Fig. 4 due 
to the smaller relative gluon contribution and the reduced number of channels available for 
the 7+jet process. Note that the relative resolved to direct production is reduced here and 
the direct process is actually dominant at y\ > for px > 10 GeV and for all y\ at pt > 100 
GeV. The antishadowing peak is higher for resolved production, shown in Fig. 14, thanks 
to the gluon contribution to resolved production. 

Finally, we show the resolved to direct ratio in Fig. 17. The direct rate alone should 
be observable at y\ > —4 for Pb+Pb, yi ~ —2.5 for Ar+Ar and y 1 ~ for O+O and 
Pt > 10 GeV. Direct production dominates over all y\ by a large factor when p? > 100 
GeV. Although the rates are lower than the dijet rates, the dominance of direct 7+jet 
production implies than the nuclear quark distribution can be cleanly studied. 

The resolved 7+jet rates are shown in Table V. These rates are only slightly larger than 
the direct rates for p T < 20 GeV but drop below the direct rates at higher p T since there is 
no large growth in the number of available production channels, as in dijet production. In 
addition, the lower effective energy plays an important role here as well. 
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D. Resolved 7+hadron production 



The leading particle px distributions of jets from 7+jet production are 

d^A^+hx = 4 I dz f e ™ d6 cm r°° dk d 3 N 7 ^ dx r 1 dx 2 
dpTd 2 b J ie min sin# cm A min k dkd 2 b A min /fe x Jx 2min X2 

x £ {F?(x, /i 2 )^(x 2 , /i 2 , b, z) + Fj(x, /i 2 )^, ix\ b, z)} 



(hi) 



x 5 



fk 



D h /k(z c ,^ 2 ) 



(16) 



The subprocess cross sections, da/di, are related to s 2 da/didu in Eq. (13) through the 
momentum-conserving delta function 5(s + t + u) and division by s 2 . 

The resolved pt distributions for hadrons are shown in Fig. 18(a). Note that the re- 
solved cross section for leading hadron production is similar to direct production, shown in 
Fig. 15(a). The same effect is seen for dijet production in Figs. 10 and 6. 

The shadowing ratios are shown in Fig. 18. The difference between the shadowing ra- 
tios for pions produced by quarks and antiquarks is rather large and reflects both gluon 
antishadowing at low px as well as the relative valence to sea contributions for quark and 
antiquark production through q(q)g — > ?(?)7- In the FGS calculations, the antiquark ratio 
reflects the flattening of the antiquark and gluon ratios at x > 0.2. Since pions produced by 
gluons come from the qq — ► jg channel alone, only a small effect is seen, primarily in the 
EMC region. Now the total pion rates follow those for quark and antiquark production of 
final-state pions instead of those of the gluon. 

Although our p^-dependent calculations have focused on the midrapidity region of \yi\ < 
1, we have shown that extending the rapidity coverage could lead to greater sensitivity to 
the small x 2 region and larger contributions from direct photoproduction, especially at low 
Pt- 

Thus 7+jet production is a good way to measure the nuclear quark distribution functions. 
Direct photoproduction is dominant at central rapidities for moderate values of pt- Final- 
state hadron production is somewhat larger for direct production so that, even if the rates 
are low, the results will be relatively clean. 



IV. SUMMARY 

There are a number of uncertainties in our results. All our calculations are at leading 
order so that there is some uncertainty in the total rate, see Refs. [9,40]. Some uncertainty 
also arises from the scale dependence, both in the parton densities and in the fragmentation 
functions. The fragmentation functions at large z c also introduce uncontrollable uncertain- 
ties in the rates. Hopefully more data will bring the parton densities in the photon, proton 
and nucleus under better control before the LHC begins operation. The data from RHIC 
also promises to bring the fragmentation functions under better control in the near future. 

While the photon flux is also an uncertainty, it can be determined experimentally. The 
hadronic interaction probability near the minimum radius depends on the matter distribution 
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in the nucleus. Our calculations use Woods-Saxon distributions with parameters fit to 
electron scattering data. This data is quite accurate. However, electron scattering is only 
sensitive to the charge distribution in the nucleus. Recent measurements indicate that 
the neutron and proton distributions differ in nuclei [41]. This uncertainty in the matter 
distribution is likely to limit the photon flux determination. 

The uncertainty in the photon flux can be reduced by calibrating it with other mea- 
surements such as vector meson production, 7 A — > V A. Studies of well known two-photon 
processes, like lepton production, can also help refine the determination of the photon flux. 
With such checks, it should be possible to understand the photon flux in pA relative to AA 
to better than 10%, good enough for a useful shadowing measurement. 
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supported in part by the Division of Nuclear Physics of the Office of High Energy and Nuclear 
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TABLES 



A 


L (mb-V 1 ) 


VS NN (GeV) 


-Ebeam (GeV) 


1L 


femax (GeV) 


£max (TeV) 


v/S 7 jv (GeV) 





160 


7000 


3500 


3730 


243 


1820 


1850 


Ar 


43 


6300 


3150 


3360 


161 


1080 


1430 


Pb 


0.42 


5500 


2750 


2930 


81 


480 


950 



TABLE I. Pertinent parameters and kinematic limits for AA collisions at the LHC [9]. We first 
give the luminosities and the NN collision kinematics, the nucleon-nucleon center of mass energies, 
\/Snn, the corresponding beam energies, -Ebeam; Lorentz factors, jl- We then present the photon 
cutoff energies in the center of mass frame, /c m ax, and in the nuclear rest frame, E m3iX , as well as 
the equivalent maximum photon-nucleon center of mass energies, y/S^N- 





Pb+Pb 


Ar+Ar 




O+O 


Pt v v j 


do j dpx 


Rate 


du j dpx 


Rate 


da j dpx 




Rate 


11.4 


2.802 xlO 5 


1.177xl0 5 


4.771xl0 3 


2.052xl0 5 


5.506x10 


2 


8.810xl0 4 


21.4 


1.389 xlO 4 


5.836 xlO 3 


2.869xl0 2 


1.234xl0 4 


3.706x10 


1 


5.930xl0 3 


31.4 


1.907xl0 3 


8.010xl0 2 


4.574X10 1 


1.967xl0 3 


6.462x10 





1.034xl0 3 


41.4 


4.181xl0 2 


1.756 xlO 2 


1.137X10 1 


4.889xl0 2 


1.734x10° 


2.774xl0 2 


51.4 


1.207xl0 2 


5.069X10 1 


3.668x10° 


1.577xl0 2 


5.983x10" 


-1 


9.573X10 1 


61.4 


4.172X10 1 


1.752X10 1 


1.403x10° 


6.033X10 1 


2.434x10" 


-1 


3.894X10 1 


71.4 


1.640X10 1 


6.890x10° 


6.061X10- 1 


2.606X10 1 


1.113x10" 


-1 


1.781X10 1 


81.4 


7.105x10° 


2.984x10° 


2.868X10" 1 


1.233X10 1 


5.557x10" 


-2 


8.891x10° 


91.4 


3.316x10° 


1.393x10° 


1.457X10" 1 


6.265x10° 


2.970x10" 


-2 


4.752x10° 


101 


1.640x10° 


6.888X10- 1 


7.832xl0~ 2 


3.368x10° 


1.676x10" 


-2 


2.682x10° 


111 


8.538X10- 1 


3.586X10" 1 


4.416xl0~ 2 


1.899x10° 


9.898x10" 


-3 


1.584x10° 


121 


4.624X10- 1 


1.942X10" 1 


2.586xl0~ 2 


1.112x10° 


6.066x10" 


-3 


9.706X10- 1 


131 


2.588xl0 _1 


1.087X10" 1 


1.564xl0~ 2 


6.725X10- 1 


3.835x10" 


-3 


6.136X10- 1 


141 


1.493X10" 1 


6.270xl0~ 2 


9.737xl0~ 3 


4.187X10" 1 


2.493x10" 


-3 


3.989X10- 1 


151 


8.838X10" 2 


3.712xl0~ 2 


6.215xl0~ 3 


2.672X10" 1 


1.660x10" 


-3 


2.656X10" 1 


161 


5.346xl0~ 2 


2.245xl0~ 2 


4.050xl0~ 3 


1.741X10" 1 


1.128x10" 


-3 


1.805xl0 _1 


171 


3.293xl0~ 2 


1.383xl0" 2 


2.687xl0~ 3 


1.155X10" 1 


7.803x10" 


-4 


1.248X10" 1 


181 


2.064xl0~ 2 


8.669xl0~ 3 


1.813xl0" 3 


7.796 xlO" 2 


5.485x10" 


-4 


8.776xl0~ 2 


191 


1.316xl0~ 2 


5.528xl0~ 3 


1.244xl0~ 3 


5.349xl0~ 2 


3.917x10" 


-4 


6.267xl0~ 2 


201 


8.507xl0" 3 


3.573xl0" 3 


8.642xl0~ 4 


3.716xl0~ 2 


2.833x10" 


-4 


4.533xl0" 2 



TABLE II. The differential cross sections, in units of nb/GeV, without shadowing and rates 
for | j/i | < 1 in a one month, 10 6 s, run for direct dijet photoproduction in peripheral AA collisions. 
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Pb+Pb 


Ar- 


-Ar 




O+O 


Pt (GeV) 


da/dpT 


Rate 


da/dpT 


Rate 


da 1 dpT 




Rate 


11.4 


5.184xl0 5 


2.177xl0 5 


1.185xl0 4 


5.096xl0 5 


1.747xl0 3 


2.795xl0 5 


21.4 


1.733 xlO 4 


7.279 xlO 3 


4.713xl0 2 


2.027xl0 4 


7.715x10 


i 


1.234xl0 4 


31.4 


1.895xl0 3 


7.959 xlO 2 


5.927X10 1 


2.549xl0 3 


1.054x10 


i 


1.686xl0 3 


41.4 


3.530xl0 2 


1.483 xlO 2 


1.247X10 1 


5.362xl0 2 


2.386x10 





3.818xl0 2 


51.4 


8.945X10 1 


3.757X10 1 


3.530x10° 


1.518xl0 2 


7.210x10" 


-i 


1.154xl0 2 


61.4 


2.767X10 1 


1.162X10 1 


1.210x10° 


5.203X10 1 


2.629x10" 


-i 


4.206X10 1 


71.4 


9.857x10° 


4.140x10° 


4.754X10- 1 


2.044X10 1 


1.094x10" 


-l 


1.750X10 1 


81.4 


3.905x10° 


1.640x10° 


2.068X10- 1 


8.892x10° 


5.031x10" 


-2 


8.050x10° 


91.4 


1.679x10° 


7.052X10- 1 


9.726xl0~ 2 


4.182x10° 


2.496x10" 


-2 


3.994x10° 


101 


7.687X10- 1 


3. 229x10^ 


4.864xl0~ 2 


2.092x10° 


1.315x10' 


-2 


2.104x10° 


111 


3.725X10- 1 


1.565X10" 1 


2.566xl0~ 2 


1.103x10° 


7.296x10- 


-3 


1.167x10° 


121 


1.884X10- 1 


7.913xl0~ 2 


1.411X10" 2 


6.067X10- 1 


4.215x10- 


Q 
O 


6.744X10- 1 


131 


9.870xl0~ 2 


4.145xl0~ 2 


8.036xl0~ 3 


3.455X10" 1 


2.519x10- 


-3 


4.030X10" 1 


141 


5.353x10 z 


2.248x10 


4.728 xlO" 3 


2.033X10" 1 


1.554x10- 


-3 


2.486x10 


lol 


2.986xl0~ 2 


1.254xl0~ 2 


Z.oOo X 1U 


l.zzy X 1U 


9.840x10- 


-4 


1.574xl0 _1 


161 


1.704xl0~ 2 


7.157xl0~ 3 


1.767xl0" 3 


7.598xl0" 2 


6.372x10" 


-4 


1.020xl0 _1 


171 


9.929xl0~ 3 


4.170xl0" 3 


1.115xl0~ 3 


4.795xl0~ 2 


4.207x10" 


-4 


6.731xl0~ 2 


181 


5.895xl0~ 3 


2.476 xlO" 3 


7.159xl0" 4 


3.078xl0~ 2 


2.827x10" 


-4 


4.523xl0~ 2 


191 


3.569xl0~ 3 


1.499xl0~ 3 


4.685xl0~ 4 


2.015xl0~ 2 


1.935x10" 


-4 


3.096xl0~ 2 


201 


2.192xl0" 3 


9.206xl0" 4 


3.110xl0" 4 


1.337xl0" 2 


1.342x10" 


-4 


2.147xl0" 2 



TABLE III. The differential cross sections, in units of nb/GeV, without shadowing and rates 
with \yi\ < 1 in a one month, 10 6 s, run for resolved dijet photoproduction in peripheral AA 
collisions. 
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Pb+Pb 




Ar- 


-Ar 






0+0 




Pt (GeV) 


da 1 dpx 




Rate 




da/dp T 




Rate 




do - / (ipy 




Rate 




11.4 


1.877x10 


2 


7.886x10 


l 


2.838x10° 


1.220x10 


l 


3.024x10" 


-l 


4.838x10 


l 


21.4 


1.435x10 


1 


6.026x10 





2.534x10" 


-l 


1.090x10 


l 


2.948x10" 


-2 


4.716x10 





31.4 


1.656x10 





1.116x10 





5.310x10" 


-2 


2.284x10 





6.630x10 


-3 


1.061x10 





41.4 


7.326x10" 


-1 


3.076x10" 


-i 


1.633x10" 


-2 


7.022x10" 


-l 


2.168x10" 


-3 


3.468x10" 


-1 


51.4 


2.546x10" 


-1 


1.069x10" 


-i 


6.264x10" 


-3 


2.694x10" 


-l 


8.804x10" 


-4 


1.408x10" 


-1 


61.4 


1.027x10" 


-1 


4.314x10" 


-2 


2.778x10" 


-3 


1.195x10" 


-l 


4.116x10' 


-4 


6.586x10" 


-2 


71.4 


4.602x10" 


-2 


1.933x10" 


-2 


1.362x10" 


-3 


5.856x10" 


-2 


2.124x10" 


-4 


3.400x10" 


-2 


81.4 


2.228x10" 


-2 


9.358x10" 


-3 


7.200x10" 


-4 


3.096x10" 


-2 


1.179x10" 


-4 


1.887x10" 


-2 


91.4 


1.146x10" 


-2 


4.814x10" 


-3 


4.034x10" 


-4 


1.735x10" 


-2 


6.930x10 


-5 


1.109x10" 


-2 


101 


6.176x10" 


-3 


2.594x10" 


-3 


2.366x10" 


-4 


1.017x10" 


-2 


4.258x10" 


-5 


6.812x10" 


-3 


111 


3.464x10" 


-3 


1.455x10" 


-3 


1.441x10" 


-4 


6.196x10" 


-3 


2.716x10" 


-5 


4.346x10" 


-3 


121 


2.006x10" 


-3 


8.426x10" 


-d 


9.052x10" 


o 


3.892x10" 


-3 


1.785x10" 


-5 


2.856x10" 


-3 


131 


1.192x10" 


-3 


5.004x10" 


-4 


5.832x10" 


-5 


2.508x10" 


-3 


1.203x10" 


-5 


1.924x10" 


-3 


141 


7.254x10 


-4 


3.046x10 


-4 


3.846x10" 


-5 


1.654x10" 


-3 


8.292x10" 


-6 


1.327x10 


-3 


lol 


4.506x10" 


-4 


1.893x10" 


-4 


2.586x10" 


-5 


1 1 10v1H" 

l.llz X 1U 


-3 


r oof? v in- 
O.OZDX 1U 


-6 


9.322x10" 


-4 


161 


2.848x10" 


-4 


1.196x10" 


-4 


1.769x10" 


-5 


7.604x10" 


-4 


4.160x10" 


-6 


6.656x10" 


-4 


171 


1.825x10" 


-4 


7.666x10" 


-5 


1.227x10" 


-5 


5.274x10" 


-4 


3.014x10" 


-6 


4.822x10" 


-4 


181 


1.186x10" 


-4 


4.982x10" 


-5 


8.620x10" 


-6 


3.706x10" 


-4 


2.210x10" 


-6 


3.536x10" 


-4 


191 


7.816x10" 


-5 


3.282x10" 


-5 


6.134x10" 


-6 


2.638x10" 


-4 


1.642x10" 


-6 


2.628x10" 


-4 


201 


5.204x10" 


-5 


2.186x10" 


-5 


4.412x10" 


-6 


1.897x10" 


-4 


1.233x10" 


-6 


1.972x10" 


-4 



TABLE IV. The differential cross sections, in units of nb/GeV, without shadowing and rates 
in the interval \yi\ < 1 in a one month, 10 6 s, run for direct 7+jet photoproduction in peripheral 
A A collisions. 
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Pb+Pb 




Ar- 


-Ar 






0+0 




Pt (GeV) 


da 1 dpT 




Rate 




da/dpT 




Rate 




da 1 dpT 




Rate 




11.4 


1.154x10 


2 


4.847x10 


l 


2.456x10° 


1.056xl0 2 


3.333x10" 


-1 


5.333x10 


l 


21.4 


5.352x10 





2.248x10 





1.391x10" 


-l 


5.981x10 





2.137x10" 


-2 


3.419x10 





31.4 


7.037x10" 


-1 


2.956x10" 


-i 


2.126x10" 


-2 


9.142x10" 


-l 


3.588x10" 


-3 


5.741x10" 


-i 


41.4 


1.497x10" 


-1 


6.287x10" 


-2 


5.123x10" 


-3 


2.203x10" 


-l 


9.352x10- 


-4 


1.496x10" 


-i 


51.4 


4.216x10" 


-2 


1.771x10" 


-2 


1.611x10" 


-3 


6.927x10" 


-2 


3.151x10- 


-4 


5.042x10" 


-2 


61.4 


1.427x10" 


-2 


5.993x10" 


-3 


6.037x10" 


-4 


2.596x10" 


-2 


1.257x10- 


-4 


2.011x10" 


-2 


71.4 


5.502x10" 


-3 


2.311x10" 


-3 


2.561x10" 


-4 


1.101x10" 


-2 


5.647x10- 


-5 


9.035x10" 


-3 


81.4 


2.338x10" 


-3 


9.820x10" 


-4 


1.192x10" 


-4 


5.126x10" 


-3 


2.776x10- 


-5 


4.442x10" 


-3 


91.4 


1.071x10" 


-3 


4.498x10" 


-4 


5.962x10" 


-5 


2.564x10" 


-3 


1.463x10- 


-5 


2.341x10" 


-3 


101 


5.198x10" 


-4 


2.183x10" 


-4 


3.155x10" 


-5 


1.357x10" 


-3 


8.144x10" 


-6 


1.303x10" 


-3 


111 


2.655x10" 


-4 


1.115x10" 


-4 


1.752x10" 


-5 


7.534x10" 


-4 


4.749x10- 


-6 


7.598x10" 


-4 


121 


1.410x10" 


-4 


5.922x10" 


-5 


1.011x10" 


o 


4.347x10" 


_4 


2.875x10- 


-6 


4.600x10" 


-4 


131 


7.736x10" 


-5 


3.249x10" 


-5 


6.025x10" 


-6 


2.591x10" 


-4 


1.795x10- 


-6 


2.872x10" 


-4 


141 


4.376x10 


_5 


1.838x10 


-5 


3.696x10" 


-6 


1.589x10" 


-4 


1.153x10- 


-g 


1.845x10" 


-4 


lol 


2.539x10" 


-5 


1.066x10" 


-5 


2.323x10" 


-6 


9.989x10" 


-5 


( .ooy X 1U 


-7 


1.214x10" 


-4 


161 


1.505x10" 


-5 


6.321x10" 


-6 


1.491x10" 


-6 


6.411x10" 


-5 


5.095x10" 


-7 


8.152x10" 


-5 


171 


9.078x10" 


-6 


3.813x10" 


-6 


9.743x10" 


-7 


4.189x10" 


-5 


3.482x10" 


-7 


5.571x10" 


-5 


181 


5.571x10" 


-6 


2.340x10" 


-6 


6.471x10" 


-7 


2.783x10" 


-5 


2.419x10- 


-7 


3.870x10" 


-5 


191 


3.478x10" 


-6 


1.461x10" 


-6 


4.369x10" 


-7 


1.879x10" 


-5 


1.707x10" 


-7 


2.731x10" 


-5 


201 


2.200x10" 


-6 


9.240x10" 


-7 


2.988x10" 


-7 


1.285x10" 


-5 


1.220x10" 


-7 


1.952x10" 


-5 



TABLE V. The differential cross sections, in units of nb/GeV, without shadowing and rates in 
the interval \y\\ < 1 in a one month, 10 6 s, run for resolved 7+jet photoproduction in peripheral 
A A collisions. 
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FIGURES 




(g) (h) 



FIG. 1. Feynman diagrams for dijet photoproduction from direct, (a) and (b), and resolved 
photons, (c)-(h). Only a sample of the resolved diagrams are shown. Crossed diagrams are not 
shown. 
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FIG. 2. The EKS98 and FGS shadowing parameterizations are compared at the scale fx = 10, 
100 and 400 GeV for (a) valence quarks, (b) sea quarks and (c) gluons. The solid curves are the 
EKS98 parameterization, the dashed, FGS. The lower values of \x give the lowest values of S\. 
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100 200 300 400 100 200 300 400 



P T (GeV) p T (GeV) 

FIG. 3. Direct jet photoproduction in peripheral collisions, (a) The pr distributions for \yi\ < 1 
are shown for A A collisions. The solid curves is the total for Pb ions while the produced quarks 
(dashed), antiquarks (dotted) and gluons (dot-dashed) are shown separately. The total production 
for Ar (dot-dot-dot-dashed) and O (dot-dash-dash-dashed) ions are also shown, (b) The fraction of 
gluon-initiated jets as a function of px for Pb+Pb (solid), Ar+Ar (dashed) and 0+0 (dot-dashed) 
interactions, (c) The EKS98 shadowing ratios for produced partons. The solid curve is the total for 
Pb ions while the ratios for produced quarks (dashed), antiquarks (dotted) and gluons (dot-dashed) 
are shown separately. The total ratios for Ar (dot-dot-dot-dashed) and O (dot-dash-dash-dashed) 
ions are also shown, (d) The same as (c) for FGS. 
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yi yi 

FIG. 4. We compare the rapidity distributions of direct and resolved dijet production without 
shadowing in peripheral collisions. The left-hand side shows the results for pt > 10 GeV for (a) 
Pb+Pb, (c) Ar+Ar and (e) O+O collisions while the right-hand side is for pt > 100 GeV for 
(b) Pb+Pb, (d) Ar+Ar and (f) O+O collisions. The solid curves are the direct results while the 
dashed curves show the resolved results. The photon is coming from the left. 
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0.8 — 



(a) Pb+Pb 






0.0 2.5 5.0 -5.0 -2.5 0.0 

yi yi 

FIG. 5. We compare shadowing ratios in direct and resolved dijet production in peripheral 
collisions. The left-hand side shows the results for pt > 10 GeV for (a) Pb+Pb, (c) Ar+Ar and 

(e) O+O collisions while the right-hand side is for px > 100 GeV for (b) Pb+Pb, (d) Ar+Ar and 

(f) O+O collisions. The solid and dashed curves give the direct ratios for the EKS98 and FGS 
parameterizations respectively. The dot-dashed and dotted curves show the resolved ratios for 
the EKS98 and FGS parameterizations respectively. The photon comes from the left. Note the 
difference in the y-axis scales here. 
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Pt ( GeV ) P T (GeV) 

— r 1 1 1 1 I 1 1 — 1 1— 




i I i I i I i I I i I i I i I i 

100 200 300 400 100 200 300 400 

P T (GeV) p T (GeV) 

FIG. 6. Direct photoproduction of leading hadrons in peripheral collisions, (a) The px dis- 
tributions for \y\\ < 1 are shown for AA collisions. The solid curve is the total for Pb+Pb 
while the produced pions (dashed), kaons (dot-dashed) and protons (dotted) are shown separately. 
The total production for Ar+Ar (dot-dot-dot-dashed) and 0+0 (dot-dash-dash-dashed) are also 
shown, (b) The fraction of gluon-initiated hadrons as a function of pr- The curves are the same 
as in (a), (c) The EKS98 shadowing ratios for produced pions. The solid curve is the total for 
Pb+Pb while the ratios for pions produced by quarks (dashed), antiquarks (dotted) and gluons 
(dot-dashed) are shown separately. The total ratios for Ar+Ar (dot-dot-dot-dashed) and 0+0 
(dot-dash-dash-dashed) are also shown, (d) The same as (c) for FGS. 
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FIG. 7. The average value of the nucleon parton momentum fraction x as a function of trans- 
verse momentum. Results are given for (a) direct and (b) resolved gluon jet production and for (c) 
direct and (d) resolved pion production by gluons. The results are given for O+O (dot-dashed), 
Ar+Ar (dashed) and Pb+Pb (solid) interactions. 
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FIG. 8. Resolved dijet photoproduction in peripheral AA collisions, (a) The Pb+Pb jet pr 
distributions with \yi\ < 1 are shown for quarks (dashed), antiquarks (dotted), gluons (dot-dashed) 
and the total (solid). We also show the total jet pr distributions in Ar+Ar (dot-dot-dot-dashed) 
and 0+0 (dash-dash-dash-dotted) collisions, (b) The relative EKS98 shadowing contributions 
from quarks (dashed), antiquarks (dotted) and gluons (dot-dashed) as well as the total (solid) are 
shown for Pb+Pb collisions. The totals are also shown for Ar+Ar (dot-dot-dot-dashed) and 0+0 
(dash-dash-dash-dotted) interactions, (c) The same as (b) for FGS. 
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FIG. 9. We present the resolved/direct dijet production ratios as a function of rapidity. In (a) 
we show the results for pt > 10 GeV while in (b) we show px > 100 GeV. The curves are Pb+Pb 
(solid), Ar+Ar (dashed) and O+O (dot-dashed). The photon comes from the left. 
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FIG. 10. Resolved leading hadrons from dijet photoproduction in peripheral collisions, (a) 
The pr distributions for \yi\ < 1 are shown for AA collisions. The Pb+Pb results are shown for 
charged pions (dashed), kaons (dot-dashed), protons (dotted) and the sum of all charged hadrons 
(solid). The charged hadron p? distributions are also shown for Ar+Ar (dot-dot-dot-dashed) and 
0+0 (dot-dash-dash-dashed) collisions, (b) The EKS98 shadowing ratios for produced pions. For 
Pb+Pb collisions, we show the ratios for pions produced by quarks (dashed), antiquarks (dotted), 
gluons (dot-dashed) and the total (solid) separately. The ratios for pions produced by all partons 
are also shown for Ar+Ar (dot-dot-dot-dashed) and 0+0 (dot-dash-dash-dashed) collisions, (c) 
The same as (b) for FGS. 
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(c) 

FIG. 11. Feynman diagrams for 7+jet photoproduction from direct, (a), and resolved photons, 
(b) and (c). Crossed diagrams are not shown. 
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FIG. 12. Direct 7+jet photoproduction in peripheral collisions, (a) The px distributions for 
I j/i I < 1 are shown for Pb+Pb (solid), Ar+Ar (dashed) and O+O (dot-dashed) collisions, (b) The 
EKS98 shadowing ratios are shown for Pb+Pb (solid), Ar+Ar (dashed) and 0+0 (dot-dashed) 
while the corresponding FGS ratios are shown for Pb+Pb (dotted), Ar+Ar (dot-dot-dot-dashed) 
and 0+0 (dot-dash-dash-dashed) collisions. The photon comes from the left. 
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FIG. 13. We compare the rapidity distributions of direct and resolved 7+jet photoproduction 
in peripheral collisions. The left-hand side shows the results for p? > 10 GeV for (a) Pb+Pb, (c) 
Ar+Ar and (e) O+O collisions while the right-hand side is for px > 100 GeV for (b) Pb+Pb, (d) 
Ar+Ar and (f) O+O collisions. The solid curves are the direct results while the dashed curves 
show the resolved results. The photon comes from the left. Note the different scales on the y-axes. 
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FIG. 14. We compare shadowing ratios in direct and resolved 7+jet production in peripheral 
collisions. The left-hand side shows the results for pt > 10 GeV for (a) Pb+Pb, (c) Ar+Ar and 

(e) O+O collisions while the right-hand side is for px > 100 GeV for (b) Pb+Pb, (d) Ar+Ar and 

(f) O+O collisions. The solid and dashed curves give the direct ratios for the EKS98 and FGS 
parameterizations respectively. The dot-dashed and dotted curves show the resolved ratios for the 
EKS98 and FGS parameterizations respectively. 
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FIG. 15. Direct leading hadrons from 7+jet photoproduction in peripheral collisions, (a) The 
Pt distributions for \yi\ < 1 are shown for AA collisions. The Pb+Pb results are shown for charged 
pions (dashed), kaons (dot-dashed), protons (dotted) and the sum of all charged hadrons (solid). 
The charged hadron px distributions are also shown for Ar+Ar (dot-dot-dot-dashed) and 0+0 
(dot-dash-dash-dashed) collisions, (b) The EKS98 shadowing ratios for produced hadrons. The 
results for pions, kaons and the charged hadron total (solid) are nearly identical. The proton result 
(dotted) is lower. The total charged hadron ratios for Ar+Ar (dot-dot-dot-dashed) and 0+0 
(dot-dash-dash-dashed) collisions are also shown, (c) The same as (b) for FGS. 
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FIG. 16. Resolved 7+jet photoproduction in peripheral AA collisions, (a) The Pb+Pb jet pr 
distributions with \yi\ < 1 are shown for quarks (dashed) , antiquarks (dot-dashed) , gluons (dotted) 
and the total (solid). We also show the total jet pt distributions in Ar+Ar (dot-dot-dot-dashed) 
and 0+0 (dash-dash-dash-dotted) collisions, (b) The relative EKS98 shadowing contributions 
from quarks (dashed), antiquarks (dotted) and gluons (dot-dashed) as well as the total (solid) are 
shown for Pb+Pb collisions. The totals are also shown for Ar+Ar (dot-dot-dot-dashed) and 0+0 
(dash-dash-dash-dotted) interactions. 
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FIG. 17. We present the resolved/direct 7+jet production ratios as a function of rapidity. The 
left-hand side shows the results for px > 10 GeV while the right-hand side is for px > 100 GeV. 
The curves are Pb+Pb (solid), Ar+Ar (dashed) and O+O (dot-dashed). The photon comes from 
the left. 
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FIG. 18. Resolved leading hadrons from 7+jet photoproduction in peripheral collisions, (a) 
The pr distributions for \y\\ < 1 are shown for AA collisions. The Pb+Pb results are shown for 
charged pions (dashed), kaons (dot-dashed), protons (dotted) and the sum of all charged hadrons 
(solid) . The charged hadron px distributions are also shown for Ar+ Ar (dot-dot-dot-dashed) and 
0+0 (dot-dash-dash-dashed) collisions, (b) The EKS98 shadowing ratios for produced pions. For 
Pb+Pb collisions, we show the ratios for pions produced by quarks (dashed), antiquarks (dotted), 
gluons (dot-dashed) and the total (solid) separately. The ratios for pions produced by all partons 
are also shown for Ar+Ar (dot-dot-dot-dashed) and 0+0 (dot-dash-dash-dashed) collisions, (c) 
The same as (b) for FGS. 
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